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Abstract: We perform the first analysis of Dark Matter scenarios in a constrained model 
with Dirac Gauginos. The model under investigation is the Constrained Minimal Dirac 
Gaugino Supersymmetric Standard model (CMDGSSM) where the Majorana mass terms 
of gauginos vanish. However, i?-symmetry is broken in the Higgs sector by an explicit 
and/or effective H^-term. This causes a mass splitting between Dirac states in the fermion 
sector and the neutralinos, which provide the dark matter candidate, become pseudo-Dirac 
states. We discuss two scenarios: the universal case with all scalar masses unified at the 
GUT scale, and the case with non-universal Higgs soft-terms. We identify different regions 
in the parameter space which fulfill all constraints from the dark matter abundance, the 
limits from SUSY and direct dark matter searches and the Higgs mass. Most of these points 
can be tested with the next generation of direct dark matter detection experiments. 
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1 Introduction 

With the discovery of a particle strongly resembling the Higgs boson we have a wealth of 
new electroweak precision observables, chief among these being the new scalar’s mass. It 
has therefore become an important endeavonr to nse this data to constrain models where 
this qnantity can be calcnlated as mnch as possible, in particnlar since there has so far 
been no clear detection of other new particles at the Large Hadron Collider (LHC). The 
most important framework where the Higgs boson mass can be calcnlated from top-down 
considerations is snpersymmetry (SUSY) and, particnlarly in light of the non-observation of 
light snperpartners, it has become increasingly important to consider more general models 
beyond the Minimal Snpersymmetric Standard Model (MSSM). 

One particnlarly interesting extension of the MSSM is to add Dirac masses for the 
ganginos[l-49]. This has a nnmber of motivations from both the top down - snch as allowing 
an (approximate or exact np to gravitational corrections) R-symmetry; simpler models of 
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supersymmetry breaking; permitting N = 2 supersymmetric subsectors of the theory at 
high energies - and bottom up: they allow increased naturalness through supersoftness 
[4, 50], contain new couplings that can enhance the Higgs mass, can weaken both LHC 
search bounds [51-53] and flavour constaints [8, 54, 55]. There is thus a growing literature 
on their phenomenology (some of which is cited above) to which we refer the reader. 

In previous work [42], (some of) the present authors introduced a constrained unified 
scenario where, from just a few parameters determined at a Grand Unification (GUT) scale 
MquT) the entire spectrum of superpartners and the Higgs boson mass could be determined. 
Preliminary scans determined the general features of the model. Here we shall extend this 
work by studying the consequences for dark matter under the assumption that the universe 
has a standard thermal history (i.e. no late-time reheating etc). Since the incontrovertible 
existence of dark matter is our best direct indication of physics beyond the Standard Model, 
no new scenario is complete without such an analysis. Interestingly, it can be used as an 
additional constraint to correlate with other searches (although given our ignorance of the 
thermal history of the universe it is more difficult to provide definitive exclusions). In 
addition we will employ a more accurate calculation of the Higgs mass (up to two loop 
order) by using the latest version of SARAH and the results of [56, 57]. 

To introduce Dirac gaugino masses we must add adjoint chiral multiplets for each gauge 
group. In going beyond the MSSM, we are still therefore left with the choice of how much 
additional matter to include, if any. One motivation is that the gauge couplings should 
still predict perturbative unification without the need to fine-tune the masses of additional 
heavy states; in [42] it was found that, once two-loop corrections are taken into account, the 
best choice is perhaps the simplest: we add an extra vector-like doublet of states with the 
same gauge quantum numbers as the Higgs doublets, and two additional pairs of vectorlike 
multiplets charged only under hypercharge with charges ±1. We are then left with the choice 
of top-down motivation for the supersymmetric and supersymmetry-breaking parameters; 
broadly we have the choice: 

1. An exact R-symmetry (up to eventual breaking by negligible gravitational effects). 
The model then resembles the MRSSM with the additional vector-like electrons. 

2. Approximate R-symmetry. 

We shall take the latter choice as in [42], for several reasons: 

• The R-symmetry must be broken in nature, and, since it is a chiral symmetry, it is 
reasonable to expect that it is broken in the Higgs sector. 

• We can then write new Higgs couplings which enhance its mass at tree level. In 
contrast, in the MRSSM the one-loop corrections due to the new higgsino states tend 
to decrease its mass. 

• We can make the choice that the new states which assure unification are instead 
vectorlike leptons. This both simplifies the structure of the model (assuring that 
the Higgs phenomenology is almost identical to the more minimal models studied in 
[12, 22, 34, 58]) and gives new predictions for lepton flavour violation experiments. 
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The question then remains as to what extent the R-symmetry is approximate, i.e. 
which R-symmetry violating terms we shall allow. We shall, as in previous work and 
e.g. [5], make the assumption that the superpotential preserves R-symmetry and the main 
source of supersymmetry-breaking terms is R-symmetric, but that additional interactions 
that solve the n/B^ problem (and couple only to the Higgs) do not, and thus we allow a 
term^. 

With the above assumptions, we can construct a complete phenomenological scenario. 
In the next section, after a review of our particle content and notation, we will describe the 
important features of the low energy spectrum based on the high-energy assumptions. In 
section 3 we describe the generic predictions for dark matter density and direct detection 
that we expect. We shall find that, since the model predicts a pseudo-Dirac dark matter 
candidate, we alleviate the tension between finding the correct dark matter density and the 
lack of signals for direct detection so far compared to a simple Majorana neutralino. In 
section 4 we present the results of numerical scans which provide quantitative predictions. 

2 The model 

Our theory, as defined in [42] (to which we refer the reader for a more complete disucssion), 
consists of the MSSM field content Q, L, E, U, D, Hu, H^ augmented by adjoint superfields 
for hypercharge, SU (2) and SU (3) denoted S, T, O, and new unification fields with repre¬ 
sentations _ 


Field 

{SU{3),SU{2))y 

Ru 

Rd 

(1) 2)_i/2 
(1) 2)i/2 

Ei^2 

Ei^2 



The superpotential of our theory is^ 

W =Y:^UiQ,Hu - Y^^DiQjHd - Y^^E^L.Hd 

+ + XsS)HaHu + 2\THdTHu 

+ (l^i? + XsRS)RuRd + 2XTRRuTRd + {iJ-Eij + 

+ XsLRiSLiRd + 2XTLRiLiTRd + XsEijSEiEj 
- Y^.RuHdEi - Y^.RdHuEi 

~ ^LFV^i ■ ^dEj - Y^pyRuHdEj. ( 2 . 2 ) 

^We shall later also consider for phenomenological reasons the possibility of a supersymmetric singlet 
mass term. 

^We use throughout this paper the conventions of Ref. [42]. Note that the definition of the electoweak 
part of the superpotential differs from that in [12] and, furthermore, from the model file implemented in 
SARAH. To use all given numbers as input for SARAH/SPheno, the following shifts must be applied 







We retain only R-symmetry preserving soft terms, in particular the Dirac gaugino masses 
miD, rn 2 D, TTT'SD and adjoint soft masses ms,mT,Tno, except for a small term, assumed 
to come from a gravitational sector in line with the discussion in the introduction. We expect 
that in explicit models other small R-symmetry breaking terms may also be generated but 
that these will not affect the phenomenology; furthermore the inclusion of a term is 
renormalisation running consistent in that it does not lead to the generation of other classes 
of soft terms^. 

At the GUT scale, we set unihed values for the gaugino masses to be rriDo and for 
scalars mo, with the exception of adjoint scalars - and, in the following, we shall consider 
scenarios with both universal and non-universal Higgs masses. The adjoint scalars and 
the Higgs may have direct couplings to the supersymmetry-breaking mediation explaining 
the non-universality; however the GUT symmetry will enforce the triplet and octet scalars 
to have a unihed mass at the GUT scale ms (which may be different from the singlet 
ms depending on the completion). Regarding the holomorphic mass terms for the adjoint 
scalars Bs, Bt, Bq (corresponding to lagrangian terms C D —^BsS"^ + h.c. and similarly 
for the triplet and octet) we take their values to be parameters that we mostly set to zero 
except where explicitly stated otherwise. Since we are considering a gravity-mediation- 
inspired scenario (with potentially mixed F- and D-terms, see e.g. [18]) the generation of 
these terms will have a different origin to the Dirac gaugino mass (unlike in typical gauge- 
mediation scenarios) and so we cannot therefore make any strong claim about what relative 
size to expect for these terms. However, it is not unreasonable for them to be small in this 
context as there is no reason to expect a Giudice-Masiero-type term to generate them from 
F-terms. 

In [42] several predictions for the low energy spectrum based on these assumptions were 
found: 

• Unihcation takes place at (1.8 ± 0.4) x 10^^ GeV. 

• We have a compressed pattern of soft masses (with deviations of a few percent upon 
varying the input parameters): 

"^1^33 • "^Q33 • 

=0.16 : 0.39 : 0.77 : 0.79 : 0.83 : 0.93 : 1.02 

®It is obvious on dimensional grounds that B-terms cannot generate A-terms or Majorana masses on 
renormalisation group running, and they do not enter in the RGBs of terms (see e.g. the generic 
expressions of the RGBs given in Ref. [58, 59]). The only other RGBs that they could enter in, therefore, 
are for other B-terms and the scalar tadpole. Now, if we imagine B^ to be a background field which 
transforms under R, if we do not include any explicitly R-symmetry-breaking terms in the superpotential 
(such as a term) which would also therefore transform under R, then B^ cannot appear in any RGB for 
an R-symmetry-preserving quantity. In our theory this is true for Bs, Bt, Bq and also the scalar tadpole - 
so their RGBs do not depend on B^ to any order. We can also see this clearly at two loops in the explicit 
expressions given in [58]. On the other hand, other B-terms could be generated if they would also violate 
R; in our theory this is the case for B - » and B = . However, other than the fact that these are of 

’ Ei Ej Ei Ej ’ 

very little phenomenological importance, their presence in the RGBs must be proportional to the couplings 
^sisijy ^SEij which we are neglecting in this analysis. 
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where the ratios are normalised with respect to the common mass at the GUT scale 
rriQ. 

• Sleptons are heavy and quasi-degenerate with the first two generations of squarks. 
This is because the Dirac gaugino masses do not enter the squark RGBs. 

• The gaugino masses are in the ratio 

miD/rriDo ■ m 2 Z)/mDo : msD/mDo = 0.22 : 0.9 : 3.5. 


In this section we shall outline the anticipated impact of these conclusions for a neu- 
tralino dark matter particle in combination with the observed Higgs mass; in the following 
sections we will perform a detailed numerical study. 


2.1 Electroweak scalar sector 


Introducting the notation 


-|- 4mi£) 

~~ T B'j' -\- 4:Tn2f) 

the mass matrix for the GP even scalars in the basis {h,H, Sji, r)j} is: 


( -|- A/iS2^ 

AhS2l3C2l3 

Ahs 

Aht 

A/IS2/3C2/3 

Ml - Ahsjf^ 

^Hs 

Ant 

^hs 

Ahs 

1 \2 

^SR '^S 2 


\ Aht 

Ant 

XsXri 

'^TR + 


where we have defined: 


(2.3) 


(2.4) 


(2.5) 


which vanishes when A 5 and At take their N = 2 values [7]. We denote non-diagonal 
elements describing the mixing of Sr and states with the light Higgs /i: 


Ahs = v[V 2 Xsfl - gYmiDC 2 i 3 \ Aht = v[-V2Xti1 + 921^200213] ( 2 . 6 ) 


where 


while 


, A 5 At 


Ahs = gYmiDVS 2 i 3 , Am = —92^20082^ 


(2.7) 


( 2 . 8 ) 


stand for the corresponding mixing with heavier Higgs, H. The minimisation conditions 
for the potential give 


1 

vs = --— 
msR 


V 

1 

V 

tSu T T^Ahs 

, VT — ^ 


n 2 

rriTR 

R 2 


(2.9) 
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where we have included tadpoles 


d/\V 


dAV 


( 2 . 10 ) 


ts 


R 




The significance of the above is that the mixing between the singlet/triplet scalars 
and the light Higgs h is proportional to the singlet/triplet expectation values. Hence by 
reducing the mixing term we both allow a larger light Higgs mass and smaller expectation 
values for the adjoint scalars. While the former is obviously desirable, the latter is also 
vitally important for electroweak precision tests because a large triplet expectation value 
would split the masses of the W and Z from their Standard Model values. For instance, 
using Ap = (4.2 ± 2.7) x 10“'^ [60-64], we require 

Ap ~ 4^ < 1 X 10“^ (95%) (2.11) 

yZ 


which is satisfied for f t ^ 4 GeV. From equations (2.6) and (2.9) we see that for large 
tan/3, small p and/or Xt we have 


0.4 > 


/ m2D 
yriiTR 


f2 TeV\ 

V ^tr J ' 


( 2 . 12 ) 


If we were to take the supersoft limit of niTR = 2 m 2 _D this would require 5 TeV Winos. 
On the other hand, since C 2/3 < 0 this implies that we may have a partial cancellation of 
the mixing terms Ahs and Aht if p or Xt/Xr are substantial and negative. By making this 
choice we both enhance the Higgs mass and reduce the shift in the p parameter. In the 
following we shall take p negative and all other quantities positive which allows for a wider 
parameter space without tuning. 



Figure 1. Left: Ap for data from figure 8. All parameter points with a suitable relic density 
and right Higgs mass fit the current experimental bounds of Ap = (4.2 ± 2.7) x 10“^ [60-64] (black 
lines). Right: Dependence of Ap on Xt- Generically, larger Xt leads to smaller values of vt and 
subsequently to smaller Ap. The fixed parameters for this plot are: mo = 2.7 TeV, moo = 1-1 TeV, 
tan/3 = 2.8, = 4* 10® GeV^, ms = 750 GeV, ms = 3 TeV, As = 0.45 and p = —625 GeV. 

In figure 1 the value of Ap is shown for the points of the random scan, as the exper¬ 
imental value could lead to heavy constraints on our calculations. However we find, that 
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most of the parameter points all of the ones with the right relic density fall well within 
the experimental value of Ap = (4.2 ± 2.7) x 10“"^ [60-64]. As discussed above, this is 
due to a partial cancellation between the contribution of the Wino mass and Xt to vt and 
subsequently to Ap. 

2.2 Sfermions and gauginos 

From [42] we recall that, at the SUSY scale: 

miD/rriDo ■ "r2D/m.D0 : w-sd/w-do =0.22 : 0.9 : 3.5, (2.13) 

^^33 • "^Q33 • • ^“eu ■ =0.16 : 0.39 : 0.77 ; 0.79 : 0.83 : 0.93 : 1.02 

(2.14) 

These ratios will have signihcant consequences for dark matter: let us suppose hrstly that 
we wish to safely avoid constraints on squarks and the gluino by placing 

rriq-^ 2 >1.5 TeV, > 0.7 TeV 

rrig >1.5 TeV 

m^± >106 GeV. (2.15) 

The loop corrections to the gluino mass are of the order of 0(100) GeV, so neglecting them 
we hud 


niDo ^400 GeV. (2.16) 

However, this would yield a very light bino which we should bear in mind. 

For the coloured squarks, they receive supersoft corrections of 

dm? ~ mg +^ 4777 , 3 ^ ^ 0.6m|,o mq > O.SniDo- (2.17) 

37r 

Thus the coloured squark constraints translate into 

772-0 > 1.7 TeV OR 722770 > 1-9 TeV. (2.18) 

or more generally 1 + 0.8*^^^ >1.7 TeV. 

2.3 Adjoint scalars 

So far we ignored the adjoint scalars. If we take the couplings A 5 , Xt to be small, then 
ms,mT barely run at all; but 

777q ~0.8l7?7| — 0.36?72 o. (2.19) 

Without introducing a Bq term, this implies mj] > 0.7772o. 

Note that we will require a large mass for the triplet adjoint scalar to avoid a large 
triplet vev. Indeed, unless we take a large value for p and Xt to have a partial can¬ 
cellation (corresponding to some tuning) we require according to equation (2.12) that 
m-T > -y/ll-b TeV X moo ] if rnno = 400 GeV then we have m -£ > 2100 GeV. 
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2.4 Neutralinos 


As in the MSSM, the natural dark matter candidate for the CMDGSSM is the neutralino 
and it shall play a central role in the following. In contrast to the MSSM, however, we have 
two additional neutralinos and so the mass matrix is larger. In the (S, B, T^, W^, H^, H^) 
basis it is given by 


/ Ms 

miD 

0 

0 



miD 

0 

0 

0 


0 

0 

My 

m2D 


0 

0 

m2D 

0 

2 ^13 

-^SR 

2 


-Tc/3 



0 



T^/3 


-^S0 


0 


We have included supersymmetric masses Ms,Mt in the above for later convenience, but 
in the pure CMDGSSM these are both zero. 

Clearly, since the ratio of the bino to wino masses is so large at low energies, in the 
CMDGSSM the Wino will be almost decoupled from a dark matter perspective. We can 
therefore consider the LSP to be either a bino, higgsino or a mixture - and indeed the 
mixing will be very important. In the limit that |/r| ;:?> mi/r or vice versa the LSP will be 
a pseudo-Dirac particle and this has a signihcant effect on the dark matter density, as we 
shall see in the following. 


3 Consequences for dark matter 

In this section we shall examine analytically the expectations for the dark matter properties 
of a neutralino in the CMDGSSM, analysing in turn the pseudo-Dirac bino, higgsino and 
then a mixed state. 


3.1 Pseudo-Dirac bino 


If ^ is sufficiently large, then the lightest neutralino will be a pseudo-Dirac bino. We can 
then estimate the relic density to be [65]: 




lOvrx/ 8.7 X 10 ^^GeV 


T-2 




m 


ID 


E 


NfYf 

Mi 

f 


(3.1) 


We shall assume that the new vector-like fermions are too heavy for the bino to annihilate 
to. Then we hnd, using Xf = 20 and g^, = 96 that 


~0.2 X 


mp 

2 TeV 


rrin 


m 


DO 


(3.2) 


From the bounds above, we require a large mp (the minimum values compatible with the 
above are mp ~ 1.4 TeV, mop ^ 1300 GeV giving a bino of ~ 290 GeV) typically implying 
a rather heavy neutralino. However, this then presents three challenges for the Higgs mass: 
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first, we require large mo to give a significant stop contribution, while we then require 
rriDO > mo for the dark matter density - which poses difficulties for the singlet and triplet 
neutral scalars mixing; in particular, this will lead to m 2 D > and we would thus need 
a very large adjoint scalar mass to prevent erasing the Higgs quartic coupling. In addition, 
such large values of rriDo, while still remaining natural, may place into doubt the reliability 
of the two-loop Higgs mass calculation on a technical level. 

If we were to consider m/ro ~ 1300 GeV then to ensure that vt is not too large we 
require ms > 3900 GeV ~ although this assumes that At is small. 

Note that in [12] it was found that slepton coannihilation could provide an acceptable 
relic density - even with rather different masses for the Dirac bino and the sleptons (unlike 
the Majorana case). However, this is not reasonable to expect here given the boundary 
conditions in the UV. 


3.2 Higgsino 

For a sufficiently large value of m^o ^ higgsino could be the dark matter candidate. For 
miD ^ Ai the eigenstates would be split via mixing with the bino by an amount [12] 
~ /i 3 ^ ) for the values of interest this is always small enough to allow co-annihilation 
during freezeout, but much too large to allow for inelastic dark matter. Hence we can 
estimate the relic density using the conventional approximation for higgsino dark matter of 
[ 66 ]: 


~ 0.1 f 


TeV 


(3.3) 


However, such a large value of ^ requiring an even larger value of mio implies that the 
gluino mass would be greater than 16 TeV; while this could be considered it is beyond the 
realm of validity of our codes. 


3.3 Mixed bino-higgsino 

This is the most natural option: that /r is small and therefore the higgsino mixes with the 
bino. Here we can at first appeal to the results of [12]: if the LSP is much lighter than 
the sfermions then it was concluded that the correct relic density is obtained for a higgsino 
fraction fh^i = -|- jiVjgp of about 0.2. In that work expressions for the mixing were 

given in the limits |/i| ;:?> mio and j^j <C miTi, but these limits are not appropriate for the 
required higgsino fraction and we therefore give here more accurate expressions. Taking all 
of the paramters in the neutralino mass matrix (2.20) to be real and positive, except for 
IX which (as discussed above) we take negative, we define the mass eigenstates in terms of 
the original eigenstates via Xi = ^ijXj where the usual rotation matrix includes phases to 
yield the masses to be real and positive, which to leading order has relevant entries given 

by" 

Nil —Ni2 1 /\/ 2 , N 21 ~ N 22 ^ i/V2 

^We shall write |t| for clarity since we are taking p to be negative; strictly speaking we should write —T- 
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Ni5 ^ 


vcjs 


2y/2{iJ?‘ - m\^) V 


miDigY — tfsV^Xs) + \g\{t0gY + \/2A5) 


Nie ^ - 


vcp 


2^/2{g? - w? 


ID> L 


miDitisgv + >/2A5) + \g\{gY — t^V^Xs) 


N 25 - 


N 2 Q ^ 


tVC /3 


2\/2(^2 _ ^ 2 ^) [ 

ivcjs 

2^/2{g?‘ — m\ 


mioigY + tgV2Xs) - \g\{t^gY - V^Xs) 


— niiDitjsgY — V2Xs) + \g\{gY + tfsV^Xs) 


•'ID! L 

The two lightest mass eigenstate values are 

,,2 


mi ,2 ^rriiD 


V 


8(^2 _ j„2^) 


‘^V2gYXs\g\c2i3 + ?ti-id(2A|' + gy) 


'Im|s2/3 


8(//2 - m2^) 


(2A| - gy) 


(3.4) 


(3.5) 


where we have dehned the masses to be positive, and the upper sign on the second line 
corresponds to the first eigenvalue. We therefore see that the typical mass splitting for 
the two lightest eigenstates is 0(10) GeV and therefore coannihilation is always important 
in this model. This has significant consequences for the dark matter phenomenology, in 
that the relic abundance can be reduced but at the same time (as we shall see) the direct 
detection cross-sections will be suppressed. 

From the above we see that the higgsino fractions of the two lightest eigenstates are 
given by 


16(/u2 — m\j^Y 

16(/u2 — 


+ miy,){gY + 2A5) -|- 2miD\g\ ^V2c2i3gY^s + S2i3{gY ~ ^A^)^ 

+ miy,){gY + 2A|) -|- 2miD\g\ ^V2c2i3gY^s — S2i3{gY ~ ^A^)^ 

(3.6) 


If we take tjs = 1 then we can write 

M2c2 \2 2 M2q2 \2 2 

f I f , I (o 

A{\g\-miDY S{\g\+miDY'’ A{\g\+miDY 8(|//| - mio)^’ 


Hence we see that as Xg increases the second eigenvalue becomes the lighter and also 
has a larger higgsino component, while for small A 5 it is the other way round. If, for 
example, we take Xg = 0.5 then a higgsino fraction of 0.2 is possible even for large \g\,miD 
if \\g\ — mini — 100 GeV. Thus we see that in the GMDGSSM the higgsino fraction is 
naturally larger than in the MSSM for a given splitting between bino and higgsino masses, 
and this will also help to naturally obtain the correct relic abundance. 

Since the sfermion spectrum is heavy, the relevant interactions for the dark matter 
density and detection are with the higgs and Z, and to a lesser extent with the W and 
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charginos (these latter are relevant for t-channel annihilation processes). Concentrating on 
the Higgs and Z portals^, the interactions can be written as 



(3.8) 


where we can write 


Chx°iXj ^ ~ ]^i4:)isi3NjQ - cpNj^) 

+ '/2{\sNii — \TNi‘i){spN+ cpNj^) + (i o j) 


(3.9) 



{Ni^Nj^ — NiQNjQ). 


(3.10) 


In the CMDGSSM the Wino is heavy and so we can neglect A^is, A^i 4 . 

3.3.1 Relic density 

The interactions via the Z portal are spin dependent and complicated. However, from [67] 
we see that for the mass range of interest to us (> 100 GeV) to match the correct relic 
density via the Z portal alone would require 


\Ni 5 Nj 5 — NieNjel ~ 0.01 0.1 


(3.11) 


although the model in that reference contained no charginos and so it is not obvious if we 
can use that result here. 

On the other hand, for a purely Higgs-mediated interaction, we can approximate the 
relic density as coming from degenerate Majorana fermions while their mass differences 
are comparable or less than the temperature at freezeout Tf ~ m^j20. Using Micromegas 
for a single Majorana fermion and then extrapolating we find for masses above the top 
production threshold 



(3.12) 


For larger masses and smaller higgsino fractions we expect this to be the most important 
process: the annihilation cross-section is roughly proportional to the higgsino fraction, while 
for the Z portal it is roughly proportional to the higgsino fraction squared. Hence we give 



2 


(3.13) 


-|- 4:x/2{gY + 2A|')c2/3(7yA5mi£)|/r| 



®By ‘portals’ refer to (as usual) the mediators exchanged in the interactions. 
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If we take tp = \ and neglect subleading corrections in (|/r| — 'm\Y))lv^ we can write 


~ 0.25 X 


(|/x| -mxDf 


miD 


V 


A| + 12 A|) V200 GeV/ 

so for example for A 5 = 0.5, = 1 we find 


0/r^ ~ 0.112 X 


Ia^I -mip V ^ 
120 GeV j 


miD \2 
200 GeV/ 


(3.14) 


(3.15) 


However, as we increase rriiD, we note that the approximation is not good for |//| — muj < 
\svsp/y/2. Later on we shall consider As = 0.5 for which we require |/i| —miD > 80 GeV, in 
which case we should evaluate the couplings numerically. Taking As = 0.5, tg = 2.7 we find 
for |/i| — mi£) = 50 GeV that Ylij=i ~ 0.07. For masses somewhat greater than 

200 GeV, however, we find that the mass difference required to enhance the higgsino mixing 
becomes small enough to potentially bring chargino coannihilation into consideration, and 
then the above approximation for the relic density will no longer be valid. 


3.3.2 Direct detection 


Since our dark matter particle is always pseudo-Dirac in nature, only the lightest mass 
eigenstate survives and may have a rather small direct detection cross-section. Indeed we 
have spin-dependent and spin-independent cross-sections [ 68 ]: 

(75/= 8 X lO'^^^cm^ , fJsD = 3 X 10“^'^cm^ . (3.16) 

In this case to satisfy the bounds we should have |c^^o^o| < 0.05 and |c^^o^o| < 0.05. In 
the case of a single neutralino with no coannihilation this would pose a problem to obtain 
the correct relic density, however in our case these bounds are generically satisfied either 
because the pseudo-Dirac bino annihilates efficiently through the higgs or Z portals, or 
because there is additional coannihilation via the Z. We can give approximate formulae for 
these couplings (for both light eigenvalues, since the lightest state depends on the choice of 
parameters): 




4(^2 _ ^ 2 ^) 

V 

4(^2 _ ^ 2 ^) 


ITT-lDigY + 2 A 5 ) -|- |/i| (2V2c2l3gY + 52/3(ffy “ 2 A|)) 
fT^ioigY + ^'^1') + \g\ (2V2c2/3g'yA5 — S2/3(fi'y ~ 2A5)) 




C2/3(5y ~ ^AI') -|- 2\/2s2/35'y A 5 


(3.17) 


^X2X2 26cw'Sry(Ai^ — I 

We see that for larger values of mi//,/i the spin-dependent scattering will become very 
small. Taking = 2.71, A 5 = 0.5, mi/) = 200 GeV again we find 


^hx°2X°2 - ~ ^ 


100 GeV 
Ia^I - miD 


, c^x§x^ - - 0-01 X 


300 GeV f 100 GeV A 


\g\ 


\g\ - miD/ 


(3.18) 
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This implies that it is very difficult to satisfy both direct detection constraints and the 
correct relic density through a pure Higgs mediated interaction involving just a pseudo- 
Dirac bino; we will have to appeal to the Z portal or resonances. 


3.4 Resonances 


We have seen that imposing that our gluino should be heavier than current simplihed 
model bounds leads to a bino heavier than 90 GeV and therefore a Standard-Model-Higgs 
resonance is not possible. Therefore the available resonances are the MSSM-like Heavy 
higgs and pseudoscalar. At tree level and with small mixing with the singlet these all have 
masses near 




2B^ 

S2/3 


(3.19) 


as in the MSSM. This is preserved at loop level even with large Xs (up to 0.7) as can be 
seen from the example spectra in [42]. The widths of the heavy Higgs and pseudoscalars 
are rather large - neglecting QCD corrections the partial width to tops is [69] 




mjj) 

Amf \ 

mjj) 


(3.20) 


which are 0(30) GeV for 500 GeV Heavy Higgses; when we consider that this is comparable 
to the kinetic energy of the neutralino at freezeout we see that a large tuning is not necessary 
to sit near the resonance. On the other hand, if we want to consider the decoupling regime 
where the niA rnh then we would require gluinos of mass rrig > 16 mA which could 
potentially take us beyond the regime of validity of the codes. 


4 Numerical results 

We have implemented the model in the spectrum generator generator SARAH [70-75] as de¬ 
scribed in detail in ref. [42] . SARAH generates the routines for a precise spectrum calculation 
in SPheno [76, 77]. The latter is used to evolve the full RGBs at the two-loop level and to 
calculate all pole-masses for supersymmetric particles and the corresponding mixing matri¬ 
ces at the one-loop level in the DR scheme. The neutral scalar Higgs masses are calculated 
at two loops using an effective potential approach [56, 57] which we cross-checked with 
the diagrammatic two-loop results including the agOet corrections for Dirac gauginos plus 
the known + ottotb + {at + «t)^ contributions from the MSSM. The spectrum is then 
passed through the SLHA+ interface [78] to micrOMEGAS 4.1 [79-83] using CalcHep [84, 85] 
model files generated by SARAH. micrOMEGAS 4.1 allows the computation of the relic dark 
matter density as well as dark matter-nucleon cross sections relevant for direct detection 
experiments. In table 1 we list the Standard Model input parameters. 

In particular we consider two different realizations: 
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^em 

127.929338 

Gv 

1.16639- lO-^GeV"^ 

as 

0.11720 

Mz 

91.18760 GeV 

mb{mb) 

4.2 GeV 

mt 

172.9 GeV 

nir 

1.777 GeV 




Table 1. Input values for the Standard Model parameters taken at Mz unless otherwise specified. 

1. universal scalar masses : all scalar masses, the sfermions as well as the Higgs 
bosons, are fixed to a common value mo at the GUT scale, which is determined by 
the requirement <71 = g 2 - As in the MSSM, the values for \fi\ and are calculated 
from the tadpole-equations in the Higgs sector. Moreover, also vt and vs are cal¬ 
culated from the tadpole-equations using numerical methods. For low values of m^ 
the pseudoscalar octet can become tachyonic due to large negative loop corrections. 
Therefore we also allow in general for a non-zero Bq term, which also shifts the mass 
of the scalar octet but has a negligible effect on the rest of the spectrum. Note, that 
although the superpotential parameter Mq violates /^-symmetry, the corresponding 
soft SUSY-breaking bilinear parameter Bq does not. 

2. non-universal Higgs masses (NUHM): the sfermions and the scalar fields of 

have a common mass mo at the GUT-scale. The parameters g,, B^ and 
m^ = m^ = m^ are defined at Mqut: whereas rrig is given at the SUSY scale. We 
choose the value for m| at the SUSY scale because when defining m| at Mgjjt we 
often encountered tachyonic states, in particular in scenarios with large A 5 and Ay- 
The tadpole equations are solved for m^^ and m|^^, vs and vt- This leads in general 
to soft masses for the Higgs doublets at the GUT scale which are different to those 
of the other scalars. Moreover, in contrast to the universal case we also define A 5 
and At at the SUSY scale rather than at the GUT scale for a better control over the 
Higgs mass. 

Before we present our results, a few comments and conventions are in order: First, 
we want to stress that all calculations assume the standard thermal history of the uni¬ 
verse. However, it might well be that this is changed due to additional states appearing at 
the various breaking stages which could lead for example to additional entropy production 
which depend on the details of the UV-completion of the model, see e.g. [86-91] for related 
considerations. Secondly, these calculations do not include higher corrections for the an¬ 
nihilation cross sections which can have quite some impact on the resulting relic density 
[92-98]. Glearly these uncertainties are much larger than the ones from observation which 
are already at the percent level. Therefore we allow below for a somewhat larger range for an 
acceptable dark matter relic density compared to the current 2a of = [0.1153,0.1221] 
preferred by experiment [99]. 

The colour coding used in the scatter plots is always as follows: 

• Red points have Higgs mass in the range 

121 GeV <mh< 129 GeV and 0.05 < < 0.20 (4.1) 
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• Black points have a Higgs mass outside of this range. 

• Green points have a relic density outside of the range but Higgs mass within the 
correct range. 

• Grey points have a neutralino-nucleon cross-section that has already been excluded 
by spin-dependent and/or spin-independent direct detection measurements by LUX 
[100] and/or XENONIOO [101, 102]. 

The points with darker shades of red, green and grey have the same properties as the 
brighter red and green as well as black points but with the difference that the associated 
spin-independent direct detection cross section is accessible for the next generation of the 
XENON experiment, XENONIT [103]. Hence, dark red points are the ones which corre¬ 
spond to the most promising scenarios because they can be tested in the near future. 

4.1 Universal boundary conditions 

We start with the more constrained setup with universal Higgs masses. First, as in previ¬ 
ous work and discussed in the introduction we assume that the only i?-symmetry breaking 
parameter in the superpotential is an explicit H^-term. In a second analysis we also allow 
a non-zero Ms- This purely phenomenological choice allows us to explore different possi¬ 
bilities to obtain the correct relic density by changing some neutralino masses and mixing 
entries without significantly changing the masses and properties of the other supersymmet¬ 
ric particles or the Higgs bosons. 

4.1.1 Casel:M 5 = 0 


mo [GeV] 
mno [GeV] 
tan (3 

400 ... 3000 

400 ... 3000 

1.1 ... 40 

A 5 

Xt 

Bo [GeV2] 

-1.5 ... 1.5 

-I/V 2 ... I/V 2 

-2-10® ... 5-10® 

Gase 1: scans with Ms = 0: 

Ms 

0 

Bs [GeV2] 

-5-10® ... 5-10® 

Gase 2 : scans with Ms 7 ^ 0: 

Ms [GeV] 

0 ... 10^ 

Bs [GeV2] 

-5-10^ ... 5-10^ 


Table 2. Ranges of the varied parameters using universal GUT-scale boundary conditions. All 
input parameters are taken at the GUT scale. 

For the chosen parameter ranges as given in table 2 we find that the tadpole conditions 
imply \fi\ > 200 GeV. Thus, this is approximately the minimal expected mass of a higgsino- 
like LSP. On the other side, we find higgsino-like DM candidates with masses up to the 
TeV range with the correct abundance as can be seen in the right plot of figure 2. In the 
left plot we show the relic density versus the LSP mass. 

The majority of points in figure 2 with |/r| > 600 GeV have a LSP with a small higgsino 
fraction, but two nearly degenerate neutralinos as lightest states, a pseudo-Dirac bino. As 
discussed in section 2, the remaining spectrum is such that no s-channel resonances with 
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m^o [GeV] /i [GeV] 

Figure 2. Results of a random scan using the parameter ranges of table 2. The left plot shows 
the relic density fth'^ versus m^o and the right one gives the higgsino component of the LSP 
l-^i 5 p + l-^i6p of Xi versus /i. The colour code is as described in the text in the context of eq. (4.1) 

sfermions in the propagator are possible. The reason is that Dirac gangino masses do 
not contribnte to the rnnning of the sfermion mass terms. Hence, in this class of models 
the sfermion spectrnm is more degenerate than in case of the MSSM and, thns, the LHC 
bonnds on hrst and second generation sqnarks imply that all of them are well above 1 TeV. 
Therefore, the t-channel contribntion to the dark matter annihilation is also snppressed, 
and the mass difference to the LSP is too great, so that co-annihilation with sfermions is 
by far not effective enongh to rednce the relic density to an acceptable level. 

The very narrow, mostly green and red strip at lower valnes for the relic density in 
the left plot of hgnre 2 (starting at {rn^o ~ 250 GeV, Qh‘^ ~ 0.007) and going to (m^o ~ 
1100 GeV, ~ 0-1)) is popnlated by points featnring an almost pnre higgsino-like Dirac 
nentralino pair: for parameter regions with small rriQ and comparatively large muQ we hnd 
several points with the |/i| < ttiid, resnlting in a higgsino-like LSP. This strip corresponds 
to the one with |A^i 5 p -|- |Vi 6 p ~ 1 in the right plot. This is analogons to the focns point 
region in the CMSSM. As in the MSSM, pnre higgsinos annihilate too effectively and only 
at rather large masses of m^o ~ |Ai| > 800 GeV is the correct relic abnndance obtained, in 
accordance with eq. (3.3). For smaller valnes of |/r| a mixing with the bino-like states is 
necessary to get the correct relic density. 

In addition we hnd several parameter points where the spectrnm is snch that a m^o ~ 
mAoj2 or m^o/ 2 , i.e. an effective annihilation via an s-channel resonance dne to a heavy 
scalar or psendoscalar Higgs boson is possible. This is the well-known Higgs fnnnel which 
is also present in the GMSSM. However, dne to the extended Higgs sector it does not only 
occnr for large bnt also for small tan f3. We hnd that in general the most likely hnal state 
is a tt pair in these scenarios. 

There is no fnrther mechanism for generating the correct range for in this version 
of the model, becanse the sfermions are relatively heavy and thns, there is neither a bnlk 
region nor is stan co-annihilation possible as already mentioned. Moreover, as we hnd a 
lower bonnd on the nentralino mass of abont 200 GeV, if all observations are to be explained 
simnltaneonsly, neither the Higgs nor the Z resonance can be present. 
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Figure 3. Relic density in the mg — mno plane using universal boundary conditions. The relevant 
parameters have been fixed to: tan/3 = 6 and Bq = —1.2 • 10®GeV^. Furthermore, we have chosen 
As = 0.18, At = 0.39 at the GUT scale which roughly corresponds to As = 0.15, At = 0.52 at the 
SUSY scale. The blue dashed lines give the Higgs mass in GeV. The light (dark) grey area have 
a too large (too low) relic density whereas the black stripes correspond to parameter regions with 
the relic density within the bounds 0.05 < < 0.20. The red (dot-) dashed lines indicate the 

rays where the sum of the two masses of the nearly degenerate, lightest neutralinos are equal to 
the mass of the (next-to) lightest pseudoscalar mass (Higgs funnel). The yellow dot-dashed lines 
indicate where the masses add up to the heavy scalar Higgs mass. In the white, upper left region, 
the pseudoscalar octet becomes tachyonic despite Bq 7 ^ 0 because of the large loop corrections due 
to the heavy Dirac Gluino. 


In figure 3 we show the relic density in the mQ-moo plane hxing the remaining param¬ 
eters as indicated in the caption to exemplify our hndings in more detail. Black regions 
have a relic density in the range 0.05 < flh'^ < 0.2 whereas the light (dark) grey have values 
below (above) this range. In the upper left white area the pseudoscalar octet becomes 
tachyonic. Apart from the obvious focus point region which extends in the area with low 
niQ and large mno, one can see the distinct Higgs funnel which consists of two separate an¬ 
nihilation strips: on the right, the (co)annihilation —?• —)> tt, i,j = 1,2 causes the 

small relic abundance, on the left the intermediate heavy scalar Higgs is the reason for the 
(slightly less pronounced) funnel; as discussed in 3, these two states are almost degenerate. 
For completeness we give also the value of the SM-like Higgs mass (dashed lines). 

Figure 4 presents the neutralino-nucleon cross sections vs. the lightest neutralino mass, 
both for spin-independent and spin-dependent measurements. While the present-day spin- 
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Figure 4. Cross section of the spin-independent (left side) cross section of the dark matter candi¬ 
date with neutrons. The parameter points and the colour coding are the same as in figure 2. The 
full blue and yellow lines show the current bounds on asi from XENONIOO and LUX, respectively. 
The projections for the LUX run 2013/14 (yellow dashed line) and XENONIT (blue dashed line) 
are shown as well. The corresponding spin-dependent cross section is given in the right plot where 
the black line gives the current upper limit on spin-dependent annihilation cross sections asD,n as 
from XENON measurements. 
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Figure 5. Random scan over the model parameter space using universal boundaries and 10 > 
Ms/TeV > 0. The colour code is as in the previous figures. The figure on the left panel shows the 
relic density vs. the mass of the LSP. The two dips in correspond to parameter points where 
TO^o = rfiz!^ and 771^/2. The right figure shows the higgsino content vs. the /r parameter. 


independent measurements are not sensitive to the scenarios under consideration, the spin- 
dependent ones already cut into the parameter space, mostly for light (i.e. m^o < 500 GeV) 
higgsino-like LSPs. We also show the current LUX [100] and XENONIOO [101, 102] bounds, 
and the projections of the exclusion potential of the 2013/2014 LUX run [104] and the 
XENON upgrade [103]. Clearly for a signihcant part of the parameter space with a mixed 
bino/higgsino LSP either a signal can be detected or otherwise the corresponding parameters 
can be excluded. However, nearly pure higgsino-like LSPs are hardly covered, i.e. they will 
be neither discovered nor excluded by these experiments. 
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Figure 6. The neutralino relic density as function of the LSP mass for a variation of Ms in the 
range Ms = [0, 5.5] TeV using universal boundary conditions. Ms = 0 corresponds to the right 
end of the plot whereas larger values lead to a suppression of the neutralino mass. The other 
parameters have been chosen as in figure 3 with fixed toq = 0.6 TeV and rriDo = 1-1 TeV. The 
narrow green band shows the measured relic density within the experimental error. The exact 
positions of the resonances are indicated by vertical lines. They correspond to the annihilation via 
X = Z,h,Ai,A 2 ,A^, /i 2 (from left to right) 

4.1.2 Case 2: non-zero Ms 

We have seen so far that in the universal case only two possibilities exist to find the correct 
relic density: either a moderate-to-large higgsino fraction of the LSPs or a Higgs funnel 
with either a heavy scalar or pseudoscalar Higgs resonance. Other possibilities are highly 
constrained because the parameters in the neutralino mass matrix also have a large impact 
on other aspects of the model. A lighter higgsino-like LSP is forbidden by the minimum 
conditions of the vacuum, and lighter bino-singlino like states are ruled out because of the 
strong relation to the Gluino mass. One can try to circumvent these interplays by allowing 
for additional parameters in the neutralino mass matrix. The simplest option is to allow 
for Ms 7 ^ 0. This term breaks i?-symmetry, and is difficult to motivate from a top-down 
approach. However, it does not introduce Majorana mass terms for gauginos nor trilinear 
soft-breaking couplings for sfermions. Thus, the most interesting, phenomenological rele¬ 
vant differences compared to the (C)MSSM are kept. In addition, the impact of Ms on the 
properties of other particles, in particular the Higgs scalars, is very moderate. 

We have already seen in the discussion of the neutralino mass matrix in section 2.4, that 
the main effect of non-vanishing Ms consists of splitting the bino/singlino Dirac state into 
two Majorana particles even without the mixing with the higgsinos. The mass of the lighter 
state is roughly m\j^/Ms if Ms m,i£). Thus, depending on the mass ratio mi£)/Ms, the 
resulting neutralino state can become, in principle, arbitrarily light. Note that this state is 
a nearly pure gauge singlet with only a very small higgsino contribution as can also be seen 
on the the right plot of figure 5. Therefore, even neutralino masses 0(10 GeV) are not in 
conflict with data. In the left plot of this figure we show the results in the — m^o plane 

A1 

using our standard colour coding. Due to the small higgsino content new mechanisms are 
needed to obtain the proper relic density: we find resonances with the Z-boson or the light 
Higgs boson are viable possibilities. The two dips in this plot correspond to ut-lsp — ^71^/2 
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Figure 7. Analogue to figure 4 for the case of non-zero Ms- 


and ttilsp — Tn^j^. However, as the first possibility requires the admixture of a higgsino, 

the corresponding dip in is not as pronounced as the corresponding Higgs-mediated 

one. In figure 6 we show the relic density as a function of m^o. Here we have taken the same 

parameters as for figure 3 together with mo = 0.6 TeV and moo = 1-1 TeV. The variation 

of m^o stems from the variation 0 < M^/TeV < 5.5. The case Ms = 0 correspond to the 

largest neutralino mass and and the LSP for this region has a sizeable higgsino admixture 

of roughly 65 % allowing an effective Xi Xj AB [i,j = 1, 2) annihilation. With increasing 

Ms, this admixture is quickly suppressed and is already below 2 % for m^o < 350 GeV 

yielding an overabundant relic density. The rest of the plot is governed by spikes of which 

each implies a different resonant annihilation channel XiA? —t X —)• AB, i.e. the peaks 

occur at m^o — mxl"^ and correspond with decreasing m^o to X = h 2 , A^, A^, A^, h, Z. 

The notch close to m^o ~ 173 GeV stems from the opening of the ti final state. We want 
•^1 

to stress that figure 6 features every single mechanism for generating the observed relic 
abundance that is possible within this model for universal scalar masses at the GUT scale. 

The neutralino-nucleon cross sections are generically lower than in the former case with¬ 
out Ms, see figure 7, because of the dominating bino/singlino nature of the DM candidate. 
As a consequence, almost none of the possible configurations are excluded by current direct 
detection experiments. However, the next generation of experiments will be able to probe 
a significant portion of the parameter space allowing for a Higgs resonance: most scenarios 
with a higgsino admixture of roughly two per mille or more (or, conversely, scenarios with 
\fi\ < 1 TeV) have cross-sections within the reach of XENONIT, cf. figure 5. 

4.2 Non-universal boundary conditions 

We shall now relax the boundary conditions at the GUT scale by allowing soft terms for the 
Higgs doublets which are not identical to those of the other scalars. This can be motivated 
by assuming an underlying GUT theory such as S'O(IO) where all matter fields come from 
three generations of a 16, but Higgs fields descend from other representations. Unification 
of the gauge groups should unify the masses of the triplet and octet adjoint scalars, but 
depending on the embedding the singlet adjoint may have a different value (as considered 
in [42]). Finally, we trade the input values of and m'jj^ for /r and i.e. the fj, and 
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mo [GeV] 

1000 . 

. 6000 

ms [GeV] 

1200 . 

. 4200 

moo [GeV] 

500 . 

. 1700 

ms [GeV] 

100 . 

. 1400 

tan/3 

1.5 . 

. 3 

[GeV] 

-1000 . 

. -150 

A| + 2A2 

0.452 . 

. 0.752 

[GeV] 

200 . 

. 1200 


Table 3. Ranges of the varied parameters in the non-universal Higgs mass scenario. All values are 
GUT-scale input, except for ms, As and At which are defined at the electroweak scale. Ms, Bs 
and Bo are assumed to be zero. 
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m^o [GeV] 

Figure 8. Relic density versus m^o using non-universal Higgs masses and the parameter 
ranges of table 3. The colour code is as in figure 2. 

are free parameters, while the soft-terms for the doublets are obtained from the tadpole 
conditions. 

The parameter ranges considered are given in table 3. Note that we have chosen 
0.45^ < + 2 At < 0.75^ at the SUSY scale as this allows us better control of the mass of 

the lighter scalar Higgs boson [34]. 

In hgure 8 we show VLh? versus m^o for a random scan over these parameter ranges. At 
hrst glance we hnd similar features as in the case of universal boundary conditions (UBC). 
However, we hnd an upper limit of m^o < 700 GeV. This can be understood as follows: the 
large values of A| -|- 2At at the SUSY scale which are required for the mass of the lighter 
Higgs to be close to 125 GeV increase during the RGE evolution to the GUT scale. This 
in turn implies that the /i-parameter has also a signihcant RGE evolution ®) resulting in 
a signihcantly lower value for n at the SUSY scale, yielding the observed upper bound on 
m^o. As a consequence the higgsino content of x? is even for dominantly bino/singlino-like 
LSPs larger than in the UBG case. This is the reason why the maximal relic density found 
is an order of magnitude smaller compared to the UBG case. In hgure 9 we give the mass 
splitting between the lightest neutralinos versus m^o (right plot) and the higgsino content 
(left plot). We hnd that the mass splitting between the lightest and second lightest 
neutralino to be smaller than 60 GeV for all parameter points and smaller than 40 GeV for 
points with the correct dark matter properties, demonstrating that co-annihilation plays a 
large role accross the whole parameter space. Finally, we hnd an lower bound on m.^o of 

®This is in contrast to the CMSSM where it is only changed by a few per-cent 
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Figure 9. Am^o = m^o — m^o for the same points as in figure 8 using the same colour code. All 
points with right relic density have Am^o < 40 GeV. Parameter points which cannot be excluded 
by XENONIT all have even smaller mass differences. All points with a small higgsino fraction 
(< 0.05) that satisfy the correct density lie very close to a resonance. 
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Figure 10. Contours of the relic density in the for mo = 4.1 TeV, tan/3 = 2.71, 

= 4.7 • 10® GeV^, ms = 370 GeV, ms = 3.4 TeV, Ag = 0.5 and At = 0.37. Black regions have 
a relic density within 0.05 < < 0.2, in dark (light) grey regions the relic density is too low 

(high). Green dashed lines give the mass of the lightest Higgs boson, i.e. the one resembling the 
SM-Higgs boson. On the red (yellow) line one has (mgf^ = 2m^o). The upper white 

corner is excluded due to a chargino LSP and the lower right one due to a tachyonic Higgs state. 
In the left figure parameters are shown at the GUT scale, in the right one at SUSY scale. 


about 200 GeV if all observations should be explained simultaneously which is similar to 
the UBC case. 

In figure 10 we show the relic density in the /i-mDo-plane, analogous to the mo-Mi/ 2 - 
plane in the CMSSM, where the remaining free parameters were fixed to the values given 
in the caption. In the light (dark) grey coloured regions the relic density is too high (low), 
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Figure 11. Spin-independent annihilation cross section usi and spin dependent annihilation cross- 
section with neutrons asD,n of the dark matter candidate for the data from figure 8 using the same 
colour code. The blue and yellow lines show the current bounds on asi from XENONIOO and 
LUX respectively. Dashed lines are projections for the LUX run 2013/14 and the currently under 
construction XENONIT. Almost all points with correct Higgs mass and relic density can therefore 
be probed with the next generation of detectors. The black line in the right figure shows the 
XENONIOO bounds on asi,n- All data points fulfil the XENONIOO bounds on neutralino-proton 
cross sections. 

whereas in the black areas we hnd a relic density of 0.05 < VLh? < 0.2. The npper left white 
corner is exclnded dne to a chargino LSP and in the bottom right one the lightest Higgs 
becomes tachyonic. In the lower black region one has co-annihilation between the higgsino- 
like and the bino/singlino like states. rriDo ~ 600 GeV corresponds to m^o ~ 173GeV 
and, thns for this and larger valnes the tt hnal state opens np which is the reason for the 
shift and broadening of the black band. Indeed the region with the correct relic density for 
—400 GeV ^ /r ^ —250 GeV can be at least partly described by eqnation (3.12). 

The yellow (red) line indicates where the LSP has half the mass of the second lightest 
(psendo-)scalar Higgs boson corresponding to the Higgs fnnnel which actnally lead to an 
nnderabnndant relic density here. The black npper left corner is a combination of co¬ 
annihilation and the effects of a not-too-off-shell Higgs state. 

To hnish, we show in hgnre 11 the nentralino-nncleon cross sections of the same pa¬ 
rameter points of hgnre 8 together with the cnrrent experimental limits set by LUX and 
XENONIOO. The resnlts are again similar to the UBC case bnt with one important dif¬ 
ference: nearly all points which are consistent with existing data, that means within the 
extended range dne to nnknown theoretical nncertainties, will be probed by XENONIT. 
This is again a conseqnence of the increased higgsino-content of the LSP. 

5 Conclusions 

We have investigated possible dark matter scenarios in the GMDGSSM. In this model, 
the pnre ganginos are Dirac states. Only dne to the mixing with the higgsinos, which are 
Majorana particles becanse of i?-symmetry breaking, the lightest nentralino is always a 
Majorana DM candidate with a small mass splitting to the next lightest state. Dne to the 
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heavy spectrum of superpartners the annihilation is generally dominated by the exchange 
of a SM-like Higgs scalar or Z-boson. 

We studied the case of fully universal scalar masses at the GUT scale and the case 
with non-unified masses for the Higgs doublets and adjoints. In the minimal version of the 
first case, with only as an ii-symmetry breaking parameter in the potential, we found 
a lower mass of the LSP with correct relic density of about 200 GeV. For this mass range, 
viable dark matter candidates are nearly degenerate bino/singlino states whereas the most 
efficient annihilation mechanism to suppress the relic density to the allowed level is via 
resonances with heavy scalars and pseudoscalars. Higgsino-like neutralinos with the correct 
relic density are only found for LSP masses of about 800 GeV as per the classic formula 
when other states are decoupled. 

If we allow for Mg to be non-vanishing, lighter neutralinos could comprise the observed 
DM of the universe. On the other hand, if we surrender the condition of complete unification 
in the soft-breaking scalar sector, light higgsinos can more easily be obtained and so the 
higgsino admixture of a bino/singlino pair is in general larger. As a consequence, co¬ 
annihilation between the bino LSP and higgsino-like particles is possible. 

For all cases considered we found that the current direct detection experiments are 
only slightly constraining. However, the next generation of experiments like XENONIT or 
the next run of LUX can probe a large part of the parameter space with higgsino-like DM 
candidates and |/i| < 1 TeV. Naturally, the bino/singlino option is more difficult to test 
experimentally. Therefore, the scenario with non-vanishing Mg and a light LSP would still 
not be addressed by these experiments. 
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